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We study Josephson junctions between superconductors connected through the helical edge states 
of a two-dimensional topological insulator in the presence of a magnetic barrier. As the equilibrium 
Andreev bound states of the junction are 47r-periodic in the superconducting phase difference, it was 
speculated that, at finite dc bias voltage, the junction exhibits a fractional Josephson effect with half 
the Josephson frequency. Using the scattering matrix formalism, we show that signatures of this 
effect can be seen in the finite- frequency current noise. Furthermore, we discuss other manifestations 
of the Majorana bound states forming at the edges of the superconductors. 
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Unlike ordinary insulators, topological insulators (TI) 
admit robust conducting states at their boundaries. 
These states display unique properties. For instance, 
a two-dimensional quantum spin-Hall insulator has he- 
lical edge states with up spins propagating in one direc- 
tion and down spins propagating in the other direction 
[H, [H . Signatures of these helical edge states have been 
revealed in transport measurements on HgTe/CdTe ^ 
and InAs/GaSb [4[ quantum well structures. 

A conventional superconductor (S) attached to such 
edge states induces topological superconductivity by the 
proximity effect. The resulting topological superconduc- 
tor has been predicted to support zero-energy Majorana 
bound state (MBS) at an interface with a topologically 
trivial region [5]. Majorana fermions have attracted a 
lot of attention because they are promising for topologi- 
cally protected quantum computation [6|. Indeed, a pair 
of spatially separated Majorana fermions form a Dirac 
fermion that could be used as a quantum bit. As the 
information is encoded non-locally, it would be quite im- 
mune to decoherence when the MBS are far away. 

When two topological superconductors are connected, 
a spectacular fractional Josephson effect has been pre- 
dicted SHI. The Majorana bound states localized 
on either side of the junction hybridize and form an An- 
dreev bound state with energy e(cp) = ±V^DA cos((/?/2), 
where cp is the phase difference between the two super- 
conductors, D is the transmission of the junction, and 
A is the superconducting gap. The energy e(cp) is 47r- 
periodic with respect to the superconducting phase dif- 
ference, and, in the absence of inelastic processes, the 
crossing between the two states at ip = 7r is protected 
by fermion parity |8]. Thus, if the phase is varied adi- 
abatically, the system should remain in the same state. 
The resulting Josephson current would then be given by 
/ ~ dcpe((p) oc sin((/?/2). As a result, under dc bias volt- 
age V, such a system has been predicted to manifest an 
ac Josephson effect at frequency uoj/2 = eV/h, that is 
half the Josephson frequency. 

The above prediction is an out-of-equilibrium prop- 
erty. Indeed, in equilibrium, the 27r-periodicity of the 



Josephson relation is restored as the ground state cor- 
responds to the energy level e(ip) = — \^DA\ cos((^/2)| 
yielding the Josephson current. To establish equilibrium 
an infinitesimal rate of inelastic scattering, e.g., due to 
the coupling with a bath is necessary. It is well known 
that, in conventional Josephson junctions under dc bias, 
inelastic scattering is induced by non-adiabatic transi- 
tions between Andreev bound states and the continuum 
of states above the gap. Through a mechanism known as 
multiple Andreev reflections (MAR) [9[, particles in the 
junction gain an energy eV at each traversal until they 
have acquired sufficient energy to escape into the contin- 
uum. In particular, this manifests itself by a dissipative 
current flowing through the junction at subgap voltages. 
Thus, one may wonder how robust is the prediction of a 
fractional Josephson effect. The aim of the present work 
is to address this question by computing the current and 
noise of an S/TI/S junction. 

The setup is shown in Fig. [T] [Ioj |. To proceed, we 



adopt the Landauer-Biittiker formalism 11 

-E3. The 

scattering states of the system are the eigenstates <£> T = 
v-) - where u± and v± are the electron 
and hole components associated with right (+) and left 
(— ) movers, respectively - of the Bogoliubov-de Gennes 
(BdG) Hamiltonian 



U = vpa z r z + M(x)cr x + A(x)e i(t) ^ )Tz r c 



(1) 



Here, v is the Fermi velocity, p is the momentum oper- 
ator, M(x) = M0(x)0(L — x) is a transverse exchange 
field within the junction, A(x) = A[0(—x) + 0(x — L)] 
and <j>(x) = (j)iO(-x) + <p r 6(x — L) are the amplitude and 
phase of the order parameter in the left and right leads 
(ip = <fii — <fi r ), and (Ji, Tj = x, y, z) are Pauli matrices 
acting on the spin (equivalent to +/— ) and Nambu (u/v) 
spaces, respectively. All energies are measured from the 
chemical potential. 

The current operator at the junction in terms of 
electron operators of right and left movers reads / = 
ev[^(0)^+(0) — ^1(0)^- (0)]. It can be expressed 
through the scattering states by using a Bogoliubov 



transformation, ip s (x) = ^ v [u su (x)% 
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FIG. 1: a) Schematic view of the S/TI/S junction, b) Dc 
current, c) real part, and d) imaginary part of the first har- 
monic of the ac current as a function of applied bias for various 
transparencies. Here I a = GVA/e with Gn = De 2 /h. 

where s = ± and v = {e,i,a} labels an incoming state 
with positive energy e, from the lead i = /, r, and of the 
type a = e,h. 

Due to multiple Andreev reflections [9[, scattering 
states correspond to a superposition of states with energy 
e + 2neV (n integer). For instance, the wave function of 
an incoming electron (e) with energy e from the left (/) 
lead can be written in the form 



4^(0, t) = 



( ^no + a 2n A n ^ 
A n 

V &2n B n J 



i{e+2neV)t 



(2a) 



and 



( C n \ 

(l2n+lCn 
&2n+lD n 



-i[e+(2n+l)eV]t 



(2b) 



to the left (x = 0) and the right (x = L) of the junction, 
respectively. Here, a n (e) = a(e + neV), where 



a(e) 



1 

A 



lei < A, 



e - sign(e)Ve 2 - A 2 , |e| > A, 



(3) 



is the Andreev reflection amplitude, and the prefactor 
J(e) = yj\ — |«(e)| 2 is set by the normalization of the 
wavefunction. 

The exchange field M leads to backscattering of right 
movers into left movers (and vice versa). This process is 
described by unitary scattering matrices 



5 fi 



t 

-r*t/t* 



and Sh 



-r* t* 
t* rt*/t 



(4) 



for electrons and holes, respectively, that are character- 
ized by a transmission probability D = \t\ 2 . In the short 
junction limit L <C v/A - where L is the distance be- 
tween the leads - and assuming M ^> A, one obtains 
D= [l + sinh^LM/V)]- 1 

The scattering matrices relate the coefficients 
A n , B ni C n , D n through the set of equations 



s h 



^n,0 + &2nA n 
tt2n+lD n 

tt2nBn 
&2n-\Cn-\ 



(5a) 
(5b) 



that can be solved numerically [l3| . 

Assuming that the incoming quasiparticles states are 
occupied according to the equilibrium distribution at 
temperature T in the leads, the average current I(t) = 
{1(f)) in the stationary regime [14[ takes the form 



I{t) = Y J Ine 



%2neVt 



(6) 



where 



I n = < DeVS n o - J de tann -^ j2 



l 2n 



A* + a-2nA- n 



(i 



l 2(m+n) 



n) a 2m) (^m+n^ 



•(J) 



In Fig. [TJ we plot the numerical results at T = for the 
dc current Io(V) as well as the real and imaginary part of 
the first harmonic h(V) (at the conventional Josephson 
frequency) . 

Let us first discuss the dissipative current Iq(V). For 
a transparent junction (D = 1), the current flowing 
through the junction is exactly half of the current of a 
perfectly transmitting conventional Josephson junction. 



This is not surprising as the BdG Hamiltonians for both 
systems are identical in the absence of backscattering - 
however, in conventional junctions, there are two copies 
of that Hamiltonian. As V — » 0, the dissipative current 
reaches the same value (2/7r)/ c , where I c = (e/2H)A is 
the critical current of the junction. 

At finite backscattering (D < 1), we find that Iq van- 
ishes in the limit V — » 0. This current suppression can 
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be understood due to the energy gap E g = (1 — \ /r D)A 
between the bound states and the continuum. In addi- 
tion, singularities appear at voltages eV = A/q, where q 
is an integer [15]. This is to be contrasted with conven- 
tional Josephson junctions, where singularities appear at 
voltages eV = 2A/q. These singularities are associated 
with the energy gap 2A between occupied and empty 
states that quasiparticles have to overcome. Namely, new 
channels for charge transfer open at the specific voltages 
eV = 2 A/q when quasiparticles can overcome the energy 
gap by performing q — 1 Andreev reflections. In S/TI/S 
junctions, the presence of the MBS reduces the energy 
gap between occupied and empty states to A. 

This manifests itself most clearly in the tunneling 
regime (D <C 1), where one sees a current onset at 
eV = A, cf. Fig. [TJ The analytic expression for the 
dc current in the tunneling regime reads iQ Un (V) = 
(eD/h)Jdeu(e)u(e - eV)[f(e - eV) - /(e)]. Here, /(e) 
is the Fermi distribution, and the normalized density of 
states on either side of the junction, v(e) = 7rA(5(e) + 
0(|e| - A)Ve 2 - A 2 /|e|, shows the contributions of MBS 
and continuum. At T = 0, the current i^ 1111 is the sum 
of two terms. The first term corresponds to the transfer 
from the continuum to the MBS for voltages eV > A, 
whereas the second term corresponds to the (conven- 
tional) transfer from continuum to continuum for volt- 
ages eV > 2 A. The suppression of the square-root sin- 
gularity at the gap edge in v(e) explains that the singular 
behavior of Iq(V) at eV = 2 A is smooth. 

For completeness, we provide the analytical expression 
for the excess current, 7 exc = Iq — (e 2 /h)DV at eV ^> A, 



7 exc = (2e/h)(D 2 A/R)(l - (D/VR) arctan y/R), where 
R= 1 — D, which varies between in the tunneling limit 
and (8/37r)/ c at perfect transmission. 

We now turn to the ac components of the current. Mul- 
tiple Andreev reflections yield ac components at multi- 
ples of the Josephson frequency uoj. In particular, con- 
centrating on the first harmonic ii, we notice that, just 
as the dc current, both its real and imaginary part show 
MAR features at eV = A/q. More strikingly, I\ van- 
ishes in the limit V — >> (except for perfectly transpar- 
ent junctions), and in the tunneling limit, I\ is strongly 
suppressed even at finite voltages. 

At first sight, these results seem somewhat paradoxical: 
Does the vanishing of I\ indicate the absence of an ac 
Josephson effect? In fact, the fractional Josephson effect 
with frequency coj/2 is absent in the formalism from the 
outset, cf. Eq. (j6|), whereas the regular Josephson effect 
with frequency uoj vanishes at small voltage. The latter is 
in stark contrast with the case of conventional Josephson 
junctions, where the stationary value of the dc current is 
recovered as V — >• [13j . 

To resolve this puzzle, we study the current-current 
correlations. In particular, we consider the (sym- 
metrized) frequency-dependent current noise, 



S(u) = J dre iujr (Si(t)Si(t^r) + 6I{t+r)6I(t)), (8) 

where SI(t) = I(t) — and the bar denotes a time 
averaging. In terms of the coefficients A n , B ni C n , D ni 
the noise can be expressed in the following form, 



f^EE / dede' J(e)J(e')l(f(e)[l-f(e')] + f(e')[l-f(e)])x (9) 

zbo; p \ 
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FIG. 2: Zero- frequency current noise for the same trans- 
parencies as in Fig. [TJ The inset shows the effective charge 
q* as a function of applied bias for D = 0.4, 0.7, and 0.9. 
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FIG. 3: Finite frequency noise of a junction with trans- 
parency D = 0.4 for various values of the applied bias. 



Here, A n = A n +5 n o/a, and all unprimed (primed) quan- 
tities are functions of e (e'). 

We first analyze the zero-frequency noise which is 
shown in Fig. [2j Taking the ratio of the noise and the dc 
current, one can define an effective charge <?* = 5/(2 Jo). 
As shown in the inset of Fig. the effective charge is 
q* ~ A/V, which has to be compared to q* ~ 2A/V in 
a conventional Josephson junction [Tg|, [l?]]- This result 
is directly related to the modified positions of the MAR 
features in the dc current discussed above. 

More interesting is the finite-frequency noise which is 
shown in Fig. [3] for various values of the applied bias. 
The most striking feature is a peak at uj = eV. The 
peak is revealed in the analytical formula for the noise 
S tun (uj) = e^ ± I^ un (V ±uj/e) coth[(eV ±uj)/2T] in the 
tunneling regime (D <C 1); it is absent at D = 1. In 
a conventional Josephson junction, the peak is absent 
at any transmission. Here it is a consequence of the 47r- 
periodicity of the Andreev bound states. This can be un- 
derstood in the following way. The two Andreev bound 
states carry a current I± = ±(e/2K)\/~D& sin(y?/2). If 
the system would remain in one of the two states indef- 
initely, one would indeed observe a fractional Josephson 
effect. However, due to the coupling between the bound 
state and the continuum in the presence of an applied 
bias, the switching probability is non-zero. As a conse- 
quence, the fractional Josephson effect disappears in the 
average current, but its signature remains visible in the 
noise spectrum. In particular, in the limit V — »• 0, where 
the switching probability is infinitesimal, the system is 
equally likely to be in either state and the current aver- 
ages to zero. At finite bias, the switching probability is 
finite, and a conventional 27r-periodic effect remains [18j |. 

To summarize, we studied the out-of-equilibrium prop- 
erties of S/TI/S junctions. The dc current shows clear 
signatures of the Majorana bound states forming in the 
junction. The expected fractional Josephson effect with 
frequency ujj/2 is absent in the average current as the ap- 



plied bias introduces relaxation processes. However, its 
signatures can be seen in the finite-frequency noise which 
displays peaks at uj = eV. We believe that our results 
are directly applicable to the S/TI/S junction studied in 
Ref. [19], where the absence of a supercurrent signature 
in the differential resistance was attributed to the lack 
of coherence in the junction. Frequency-dependent noise 
experiments [io| could reveal whether coherent effects do 
exist in this sample. 
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